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1. Introduction 

We present the outcome of our investigation of an exactly solvable model glass 
that shows all of the basic features of much more complicated real glasses, such 
as aging j^, ^, diverging relaxation time to equilibrium (Vogel-Fulcher-Tammann- 
Hesse (VFTH) like |^), configurational entropy satisfying the Adam-Gibbs relation 
, Kauzmann transition Q and violation of Fluctuation-Dissipation Theorem ^ . 
In studying the model we make use of a particular parallel Monte Carlo dynamics 
that retains the fundamental collective nature of the glassy dynamics and that can 
be carried out analytically on our model. We will see how this dynamics undergoes a 
huge slowing down as the system is cooled down and which kind of aging dynamics 
the system sets out. We can implement such dynamics even below the Kauzmann 
temperature, thus getting information in a regime where few analytic results are 
known. 

The main motivations for this work are twofold. Firstly we want to get more 
insight in the glassy dynamics, in its various aspects, exploiting the analytical 
solvability of our model. Indeed, every detected feature of the glassy behaviour can 
be connected in a direct correspondence with given elements of the model, thanks to 
its simplicity. We can even switch on and off certain properties or certain dynamic 
behaviours, tuning the model parameters or implementing the facilitated dynamics in 
alternative ways. Furthermore, the configurational entropy is exactly computable (see 
section ^) as a function of the dynamic variables of the model. 

The second goal is to check the generality of the concept of effective temperature, 
very often discussed in literature in many different approaches (see for instance in 
Hi 01 §7 ^ [l^i |l^i [l^l ) and to see whether the possibility exists of inserting such 
a parameter in the construction of a consistent out of equilibrium thermodynamic 
theory. Even though the physics of our model is simple, we shall find general aspects 
of the results by formulating them in the thermodynamic language (see section ^) . 

In glasses, the exponential divergence of time-scales (as opposed to the algebraic 
divergence in standard continuous phase transitions), might induce an asymptotic 
decoupling of the time-decades [Q. Exploiting this, we build the dynamics of our 
model on the reasonable assumption that, in a glassy system that has aged a long 
time t, all processes with equilibration time much less than t (the f3 processes) are in 
equilibrium while those evolving on time scales much larger than t are still quenched, 
leaving the processes with time scale of order t (i.e., the a processes), as the only 
interesting ones. To model this, we will introduce by hand two different kind of 
variables evolving on two very different (decoupled) time-scales. 

Such an asymptotic decoupling of time scales, that is the input for the present 
model and for a set of other models of the same class 0, |l8), could be the 

basis for a generalization of equilibrium thermodynamics to systems out of equilibrium 

In section |^ we introduce the model and in section ^ the facilitated dynamics that 
we apply to it. The dynamical behaviour of the one-time observables will be shown 
in section ^. In section |^ an out of equilibrium formulation of thermodynamics is 
proposed. From the study of the two-time variables dynamics (in section ^) we show 
that, quite reasonably since we are not at equilibrium, the Fluctuation-Dissipation 
Theorem is not valid anymore. We will see in section ^ that from the Fluctuation- 
Dissipation Ratio we can derive an independent definition of effective temperature. 
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2. The model and its static properties 

2.1. Hamiltonian and constraint on the configuration .space 

The model displaying a fragile glass behaviour was introduced in pO[ | and widely 
studied in It is described by the following local Hamiltonian: 

Af TV JV N 

n[{x,},{S,}] = -KJ2^^ - HJ2x,~ jJ2xA~ LY,S^ , (1) 

i—1 i—1 i—1 i—1 

where N is the size of the system and {xi} and {Si} are continuous variables, the last 
satisfying a spherical constraint: Sf = N . We will call them respectively harmonic 
oscillators and spherical spins. K is the Hooke elastic constant, H is an external field 
acting on the harmonic oscillators, J is the coupling constant between {xi} and {Si} on 
the same site i, and L is the external field acting on the spherical spins. A separation 
of time scales is introduced by hand: the spins represent the fast modes and the 
harmonic oscillators the slow ones. We assume that the {Si} relax to equilibrium on 
a time scale much shorter than the one of the harmonic oscillators. From the point of 
view of the motion of the {xi}, the spins are just a noise. To describe the long time 
regime of the {xi}, we can average over this noise by performing the computation 
of the {Si} partition function, obtaining an effective Hamiltonian depending only on 
the {xi}, that determines the dynamics of these variables . Using the saddle point 
approximation for large N we find: 



Zs{{x,})^ / iT\dsAexp{-(3H[{x,}AS.}]} <5 V Sf - iV g-'^^-"-'" (2) 





where we introduce the effective Hamiltonian 'Hcff{{xi}) = — T log Zs{{xi}), that is 
the free energy for a given configuration of {xi}, as 

K TN w -\- — 

Kff({x.}) = —m2N - HmiN - wN + — log . (3) 

We defined the short-hands 



1 ^ _ 1 ^ 2 ... T2 



mi = — 2_^a;j , "^2 = 2^ a;^ , w = J^m2 + 2JLmi + + — . (4) 



=1 1=1 



The effective Hamiltonian, equation (||), can also be written in terms of the 
internal energy U{{xi}) and of the entropy Scp({xi}) of the equilibrium processes 
(i.e. the spins): 

n,si{x^}) = U{{x,}) - TS,pi{x,}) ; (5) 

Ui{x,})^^m,2N-Hm,N-wN + ^ , ^.^({^J) ^ ^ _ | fog !^i±Z^ . (6) 

The model is also characterized by a constraint on the phase space, introduced to 
avoid the existence of the single global minimum, and implementing a large degeneracy 
of the allowable lowest states. The constraint is taken on the {xi}, thus concerning 
the long time regime. It reads: 

m2- ml> mo , (7) 

where mo is a model parameter. It is a fixed, but arbitrary, strictly positive constant. 
This constraint, applied to the harmonic oscillator dynamics, is a way to reproduce 
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the behavior of good glass formers, i.e., substances for which nucleation of the crystal 
phase is especially unlikely even at very slow cooling rates (e.g., network formers 
B2O3 and Si02, molecular organics such as glycerol and atactic polystyrene, and 
different multicomponent liquid mixtures). These are substances for which there are 
non-crystalline packing modes for the particles composing them, that have intrinsically 
low energy. The amorphous configurations are thus favored. In general the crystal 
state still exists, at lower energy, but the probability of nucleating a crystal instead of 
a glass is negligible. In specific cases (binary solutions) the glassy state can even be 
lower in energy than the crystalline one and is thermodynamically stable with respect 



to any crystal configuration |21 



In the next section we will impose a Monte Carlo dynamics |15|, ^ satisfying this 
constraint and coupling the otherwise non- interacting {xi} in a dynamic way. We will 
study such dynamics analytically. 

To shorten the notation we define the modified "spring constant" K and "external 
field" H: 

k = K -, H = H+ (8) 

w + T/2' w + T/2 ^ ' 

We stress that K and H are actually functions of the {xi} themselves (through irii 
and 7712, occurring in w). 

2.2. Statics at heat-bath temperature T 

Before describing the facilitated dynamics employed we sketch very briefly the static 
of the model. The partition function of the whole system at equilibrium is: 



ZiT) = J VxVS e^p[-(3'Hi{x,},{S,})]d Ij2 



dmidm2 exp < —(3N 




K rr /w + T/2\ T, , , 

■ma - Htui - w + - log - 77 (1 + logl'^ 



2 2 ^\ T 2 



The new object (with respect to equation (^) appearing in the exponent of the above 
expression is the configurational entropy 

X = y [1 + log(m2 - m?)] . (10) 

The saddle point equations are found minimizing the expression between square 
brackets in ^ with respect to mi and m2. This yields 

g(mi,m2) _ — —2 T 

^=mi , m2 - m2 ^ - . (11) 

A(mi,m2) K[mi,m2) 

The form of the solutions mi(T), m2(T) is quite complicated because each of these 
equations is actually a fourth order equation, but they can be explicitly computed. 

3. Facilitated Monte Carlo dynamics 

We assume as dynamics a parallel Monte Carlo (MC) dynamics for the harmonic 
oscillators. This kind of analytic Monte Carlo approach was first introduced in |^^ , 
and later applied to the exactly solvable harmonic oscillator model jl^ (which is 
just our model after setting J = L = 0) and to the spherical spin model [|6[ |l^ 
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(which is the present model after setting H — K — and considering the {xi} 
as quenched random variables). The dynamical model thus obtained with a very 
simple Hamiltonian and a contrived dynamics has the benefit of being even solvable 
analytically, which yields a much deeper insight into its properties than numerical 
simulations. Moreover, in the long-time domain the dynamics looks quite reasonable 
in regard to what one might expect of any system with a VFTH-law in its approach 
to equilibrium. 

In a Monte Carlo step a random updating of the variables is performed {xi — > a; ■ = 
-\-i'i/^fN) where the {r^} have a Gaussian distribution with zero mean and variance 
A^. We define x = TL{{x'j})—TL{{xi]) as the energy difference between the new and the 
old state. If a; > the move is accepted with a probability W{fix) = exp(— else it 
is always accepted {W{(3x) = 1). The updating is made in parallel. It is the parallel 
nature of the updating that allows the collective behavior leading to exponentially 
divergent time scales in models with no interactions between particles such us ours. A 
sequential updating would not produce any glassy effect. This dynamics may induce 
glassy behavior in situations where ordinary Glauber dynamics [p3| would not. In 
our model the parallel dynamics mimics the presence of interactions between atoms 
in realistic glasses, where a large internal cooperativeness occur. In this respect the 
constraint is not essential. A parallel MC dynamics that does not take into account the 
constraint on the {xi} still maintains glassy behaviour (see for the study of such 
dynamical version of our model). Here, however, we will only look at the dynamics 
with a built-in constraint. For different examples of dynamics implying non trivial 
collective behavior the reader can look, for instance, at the n-spin facilitated kinetic 



Ising model 124, 25 or at the kinetic lattice-gas model [El, 26 



In this section we will show the basic steps leading to dynamical equations. 
For further details refer to jl^, ^ |T^. In a Monte Carlo step the quantities 
Xi = Nmi and J2i ^1 = A^'7i2 are updated. We denote their change by yi and 7/2, 
respectively. Their distribution function is, for given values of mi and m2, 

piy,,y,\m,,m,) ^ / R 7=f ^"'''^^'^'^ f E " E - J/i) (E E ^ 

i \ i i / \ i i / 

1 / yl {y2-A^-2yim,)^ \ 

4^AVm2-m^ \ SA^m^ - ml) ) ' ^ ' 

Neglecting the variations of mi and m2 of order A'^/N we can express the energy 
difference as 

K ~ 

x^ — y2~Hyi, (13) 

In terms of x and y ~ yi the distribution function can be formally written as the 
product of two Gaussian distributions: 

P{yi,y2\mi,m2)dyidy2 = dx p{x\mi,m2) dy p{y\x,mi,m2) (14) 

- cxp fJ^^Dl] cxp 

where 

x = A^K/2, A^ ^ A'^K^{m2-mj) + A'^K^ (mi- H/K^ , (15) 
y(,) = mi-H/K ^ ^x-x ^ ^ A\m2-m,l) ^ ^^^^ 



m2 



(mi - H/K^ ^ m2-ml+ (mi - H / ' 
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To represent a fragile glass the dynamics that we apply to the system is a 
generalization of the analytic treatment of Monte Carlo dynamics introduced in . 
Also in this generalized case the dynamical model with a contrived dynamics can be 
analytically solved. 

We let A^, the variance of the random updating {vi}, depend on the distance from 
the constraint, i.e. on the whole {xi} configuration before the Monte Carlo update: 



A^{t) = 8[m2{t)-ml{t)] 



B 



1 f 



(17) 



_m2{t) - m{{t) - mo_ 

where B, mo and 7 are constants. The exponent 7 also enter the VFTH relaxation 
law, as we will see later on. In literature it is usually set equal to 1, and an argument 
for this choice was given by Adam and GibbsQ. An exact explanation for it was 
provided by Kirkpatrick, Thirumalai and Wolynes in and a further quantitative 
analysis is also reported in ||2^. However, their studies do not exclude exponents 
7 > 1, always compatible with data, merely affecting the width of the fitting interval. 



Analytic approaches |29, ^ give 7 = 2 in three dimensions. Here we consider 7 as a 
model parameter, which can be chosen below, equal to, or above unity, and investigate 
aspects of this standard picture. 

The VFTH law that we obtain is a direct consequence of the special choice ( p7[ ) 
for the MC update. In the harmonic oscillator model and in the spherical spin model 
studied in pO| , 15, the dynamics was performed within this approach, but at 
fixed A. Both cases showed a relaxation time diverging at low temperature with 
an Arrhenius law, typical of strong glasses. The same is found by setting niQ = 
and 7 = 1 in the present model but here we want, instead, to develop a model 
representing a fragile glass with a Kauzmann transition at a finite temperature. 

The question whether detailed balance is satisfied is also non-trivial in our model. 
Indeed, it happens to be satisfied for this kind of dynamics only for large iV. For exact 
detailed balance we should have 

p(x|mi, 7712) exp(-/3x) = p{-x\mi,m2) (18) 

but now, when we perform the inverse move {x'^} {xi}, the probability distribution 
also depends on the {rj} through A^ as defined in equation (|l^). It can be verified 
that the violation is of order 1/A^. 

The Monte Carlo equations for the dynamics of mi and m2 derived from this 
construction read 



rhi^ J dyidy2W{Px)yip{yi,y2\mi,m2)=^ J dxW{Px)y{x)p{x\mi,m2) , (19) 
1^2= J dyidy2W{l3x)y2piyi,y2\mi,m2)^-^J dxW{l3x){x + Hy{x))p{x\mi,m2) (20) 



4. Single-time dynamical observables 

The dynamics of the system can be expressed in terms of two combinations of mi and 
m2. The first one, defined as 

fii = -^~mi. (21) 
K 

represents the distance from the instantaneous equilibrium state. By instantaneous 
equilibrium state we mean that H and K depend on the values of mi and m2 at a 
given time t. For t 00, at the true equilibrium, one has /Lti = 0. 
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The second dynamical variable is defined as the distance from the constraint : 

^2 = 'm2-m\ - Too. (22) 

When /i2 = the constraint is reached. This will happen if the temperature is low 
enough [T < To) and the time large enough. Tq is the highest temperature at which 
the constraint is asymptotically (t oo) reached by the system, it is identified with 
the Kauzmann temperature p9| . 

The nearer the system goes to the constraint (i.e. the smaller the value of 
TO2 — ml — Too), the larger the variance A^, implying almost always a refusal of the 
proposed updating. In this way, in the neighborhood of the constraint, the dynamics 
is very slow and goes on through very seldom but very large moves, which can be 
interpreted as activated processes. When the constraint is reached the variance 
becomes infinite and the system dynamics gets stuck. The system does not evolve 
anymore towards equilibrium but it is blocked in one single ergodic component of the 
configuration space. 

In terms of ^ii and /i2, the equations of motion (p^) , (|20|) become 

ill = — JQ J dx W{(3x) X p(a;|TOi, TO2)— (1 + QD) J dx W{Px) y{x) p(a;|mi, 7712), (23) 
lJ,2 — — ( dx W{(3x) xp(a;|mi, 777,2) + 2/ii ( dx W{f3x) y{x) p(a;|mi,TO2), (24) 



K 

where D and Q are given by 



D = HJ + KL = HJ + kL , Q = . ^ ,^ 2 ■ (25) 



K^w{w + T/2f 

Above To ordinary equilibrium will be achieved without reaching the constraint. 
The temperature is, then, too high for the system to notice that there is a constraint 
at all on the configurations (we are speaking about the asymptotic time regime), and 
this implies 

lim ii2{t) = Jl2{T) = - TOO > , (26) 

(T) 



K^[T)= lim K{mi{t),m2{t);T)^K{mi{T),m2{T)). (27) 

t — ^CJO 



where 



mi and TO2 are the solutions of the static self-consistent equations (|1 1[) , if T > Tq. 
For T < To the second equation in (|l^) should be replaced by 

7772 — — Too VT < To . (28) 

Below To the system goes to configurations arbitrarily close to the constraint, and 
then stay there arbitrarily long. By definition of Tq, we can write 

mo = - ^° . (29) 



The equations of motion |23| and ^ can be solved in the long time regime, for 
fixed parameters (aging setup). We notice that the value of the VFTH exponent 7 
discriminates between different dynamic regimes if7>l,7 = lorO<7<l JTot 
(the situation 7 = 1 remains model dependent even in the long time limit). We find, 
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to the leading orders of approximation for large times, the following behavior for 

li2{t) ^ ■ (30) 

[log(tAo)+clog(log(Vio))]'/^ 

For T > Tq the parameters c and to ai'e 

_1 , _ V^jl + Q^D) 

'"2' *°-87(l + P.+Q.i?) ■ 

to is of 0(1) for a large range of parameters values and is given by equation ( p5| ) 
computed at mi = mi(T) and rn2 = t^2(T). The other function of temperature 
appearing in the above expression, P^, is defined as the infinite time limit of 

p = J^(m2 - mf) 

2Kw{w + T/2f ' 

That means P computed at mi = mi(T) and TO2 = rn2{T). The solution ( |30| ) is valid 
in the aging regime, when to t Teq(T). 

Below To the qualitative behavior of ^i2{t) is the same, but T is never reached. 
Concerning the solution (^) the only difference is in the values 

(33) 



2 + 7 _ {2K^{T)mo-T 



The dynamical behavior of /ii depends not only on the temperature (above or 
below To) but also on 7 being greater, equal to or less than one. With respect to the 
relative weight of /ii and fi2 we can identify different regimes ||l9|| . In this presentation 
we just state that for the regime with T > Tq and for the one with T < Tq and 
7 > 1, /^i(t) <C fJ'2{t) and a unique effective thermodynamic parameter can be properly 
defined in various independent ways (see next section). 

What we said up to now concerns the aging regime, but, when time grows on 
even larger time scales, finally approaching equilibrium in the temperature regime 
T > Tqi the equations of motion for any one-time observable o{t) (magnetization, 
energy, distance from the constraint, etc.) take the form 

oit) . (34) 

From the study of the dynamics for very large times ^ we get a characteristic 
relaxation time to equilibrium that depends on temperature following a generalized 
VFTH law: 

A 



Teq - exp [yZTj^ J ' (^^) 

When t Teq(T) exp [A/{T - Tq)]'^ the "distance" ^2 becomes, 

B 



\T-To ) 



1/7' 



ozT-Ta . (36) 



The parameter To in the VFTH law is identified with the Kauzmann temperature, 
i.e. the temperature such that 1{To) = To is the minimum of the configurational 
entropy [and for any T < Tq remains 1{T) = Iq]. Moreover, the specific heat 
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displays a discontinuity at Tq: at that temperature the model undergoes a real 
thermodynamic phase transition. The Adam-Gibbs relation between relaxation time 
and configurational entropy density (equation ([l^)) is also achieved, in the formp9|] 



5. Effective temperature and out of equilibrium thermodynamics 

In this section we introduce effective parameters in order to rephrase the dynamics 
of the system out of equilibrium into a thermodynamic description (for a complete 
derivation see [|o)). 

In ijl^ we got, through different methods, various expressions for the effective 
temperature Tg as function of the interaction parameters of the model and of the time 
evolution of its observables. All of them were coinciding in the regime for T > Tq and 
7 > 1. 

We want to shortly recall one particular derivation of Tg- Knowing the solution 
of the dynamics at a given time t a quasi-static approach can be followed by 
computing the partition function Zg of all the macroscopically equivalent states at 
the time t. In order to generalize the equilibrium thermodynamics we assume an 
effective temperature and an effective field Hg^ and substitute the Boltzmann- 
Gibbs equilibrium measure by exp{—Tlcsi{xi},T^He)/Te), where HcB is given in (||) 
and the true external field H has been substituted by the effective field He- As we 
get the expression of the "thermodynamic" potential Fe = —Te log Ze as a function 
of macroscopic variables mi_2 and effective parameters, we can determine Tg and 
He minimizing F^ with respect to mi and m2 and evaluating the resulting analytic 
expressions at mi^2 — ^711,2 (^)- 

The partition function of the macroscopically equivalent states is: 

5{Nmi-^x,) 5{Nm2-^x1) . (38) 

i i 

From this we build the effective thermodynamic potential as a function of Te and 
He, besides of T and H , where the effective parameters depend on time through the 
time dependent values of mi and m2, solutions of the dynamics: 

Fe{t) = U- TS,p - Te{t)I +[H- HS)]Nmi{t), (39) 

with 

Te{t) =i^(mi(i),m2(t))[mo + M2(t)], (40) 

He{t) ^ H - k {m^{t),m2{t)) ^ii{t) . (41) 

Te and H^ are actually a way of describing the evolution in time of the system out of 
equilibrium. U is the internal energy of the whole system, Sep is the entropy of the fast 
or equilibrium processes (the spherical spins) (see equation (0)) while X is the entropy 
of the slow, "configurational", processes (the harmonic oscillators, see equation ([lO|)). 
The last term of Fe replaces the —HNmi occurring in U by —HeNrai. U, Sep and 
T are 'state' functions, in the sense that they depend on the state described by T, 
Te, H and, if needed, He- In the case where only one relevant effective parameter Te 
remains, these functions do not depend on the path along which its value has been 
reached. 



Vx exp 



-—'Hosi{x^},T,He) 
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As we already mentioned, for our VFTH relaxing model at T > Tq, and at T < To 
with 7 > 1, the effective temperature alone is enough for a complete thermodynamic 
description of the dominant physical phenomena (H,, = H). 

6. Two-time variables: breaking of time-translation invariance and the 
fluctuation-dissipation relation 

In this section we compute the correlation and response functions which, unlike the 
energy and the quantities mi{t) and m2{t), depend in a non-trivial way on two times, 
when the system is out of equilibrium, thus showing directly the loss of time translation 
invariance with respect to the case at equilibrium. The aim of computing such 
quantities is also to build a Fluctuation Dissipation relation and look at the meaning 
of the Fluctuation-Dissipation Ratio (FDR), dt>C{t,t')/G{t,t'), far from equilibrium, 
and to compare it with the effective temperature derived in other ways, e.g., as in 

The correlation functions between the thermodynamic fluctuation of a quantity 
ma{t) at time t and that of a quantity mi,{t') at a different time t' are defined as 

Cab{t,t') = N{5ma{t)5mi,{t')), a, 6 =1,2 (42) 

where (....) is the average over the dynamic processes, i.e., the harmonic oscillators. 

The response of an observable at time t to a perturbation in a conjugate field 
H}j at some previous time t' takes the form 

Ga,{t.t')^^-^^, a, 6 =1,2 (43) 

In our model Hi ^ H and H2~—K/2. 

In order to be concise, in the following we will only give results without derivation; 
moreover, we will concentrate on the FDR for fiuctuations of Xi — Nnii, coupled 
to the external field H, i.e., we only consider a = b = 1. For the complete derivation 
refer to [|l§. 

Knowing the evolution of the two-time observables we can generalize the 
Fluctuation Dissipation Theorem defining another effective temperature, Tf^, by 
means of the ratio between the derivative with respect to the initial time (also called 
"waiting" time) t' of the correlation function Cn and the response function Gn: 

T^^Dit t') = f44) 

Dynamics varies strongly if T is above or below the Kauzmann temperature and 
this difference produces different equations for the leading terms in the correlation 
and response functions. Therefore, we present results for the two cases separately. 

6.1. High temperature case: T > Tq, V7 

First of all we define the time evolution function for the considered time-scale sector 
as 

h{T) = exp(~ 1^ f{t)dt] (45) 



and the acceptance rate of the MC dynamics 
_ _ exD (-D T„ 
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Tg is, in the above expression, just an abbreviation for K {rriQ + 1.12) , as given in 
equation (|40|). We will look at its relation with T^^{t, t'). The function / in equation 
( p5| ) is, in this regime, 

2DQP DP{1 



f = -iT Ul + QD)T - 



1 + QD- 



QD) 

1 + QD j{l + P + QD) 



with D, Q and P given respectively by ( p5| ) and (|3 
The correlation function conies out to be 



Cnit,t') = C,^it',t')^ + 0{^il+''T) 
h{t) 



with 



Cn{t',t') 



mo + H2{t') 



l + Q{t')D ■ 

Following the approach of |lO| we also derive the response function: 



Gii(i,t') 



Gn(t',t')^ + 0(/.^+^T) 
h[t) 



with 



Gii{t',t') = ~ f3 J dyidy2W'{(3x)ylp{yi,y2\mi,m2) 

= - P j dxW'{(3x)/:^yp{x\mi,m2)+0{tilT) 
AT{t')T{t') 2T{t') 



Kit') 



Kit') 



OimT) 



Eventually we get 



Tr{t,t') 



(47) 

(48) 
(49) 

(50) 
(51) 

(52) 
(53) 



r,(0[l + O(M2(0") + O(/.2(i')')] • 
For 7 > 1 this is equal to Tg in the aging regime i^2{i'V goes to zero faster 
than /i2(i') ~ T'e(i') — T), while as soon as the VFTH exponent equals one, the 
correspondence breaks down and only the asymptotic limits of the two effective 
temperature will be equal to each other (and to the heat bath temperature). 

6.2. Low temperature case: T < Tq, 7 > 1 

Our approach also allows us to study the regime below the Kauzmann temperature 
In this last case, though, we have qualitatively different behaviors depending on the 
value of 7, i.e., on the relative weight of /^i and /X2. We describe here the case 7 > 1, 
where /xi <C ^2 fUl • 

We find the solutions of the equations for the two-time correlation functions with 
the following expressions for the function /: 



/ = -4Tr(l + QD)-8T 



QDP 



(l-3r + 2r2) + 0(Ai2T) 



(54) 



where r is an abbreviation for the normalized difference between the effective 
temperature (40) and the heat bath temperature 

-T 



2Te-T 



(55) 
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The two-time correlation function turns out to be: 



K + M2(0 + o(m2(0)]t 



h{t') 



(56) 



h{t) 



For the response function, we have 



Gn{t,t') 



4T{t')Tit') _ 2T{t'){l-2r{t'))^ 
k{t') k{t') 




(57) 



It follows that 



Tf ^(i, t') ^ Te(t') [1 + (P2{ty) + o{^,,{t'f+-<)] . 



(58) 



In this case 0(/^2) always smaller than 0(/i2): in the long time regime Tf^{t) 
coincides with Tf.{t). 

With this outcome we demonstrate that it is possible to identify an effective 
temperature that, coupled to the configurational entropy, is able to map the dynamics 
of a system out of equilibrium into a thermodynamic frame. In our model, either 
above and below the Kauzmann transition, such a construction seems to be well 
founded, provided that the VFTH exponent 7 is bigger than one. In this case the 
time evolution is so slow that a quasi-equilibrium occur at a time dependent effective 
temperature. For 7 below one, instead, the time evolution is too fast to allow for such 
a straightforward translation into thermodynamics. 

7. Conclusions 

In the present work we have been studying a particular model glass, that has all the 
basic attributes of a real glass and the dynamics of which can be analytically solved. 

An important assumption for our study has been the decoupling of time-scales 
of the processes taking place in the glassy dynamics. The decoupling of time scales is 
also fundamental for a generalization of equilibrium thermodynamics to systems far 
from equilibrium. 

From the temporal behavior of the slowly varying observables in the aging regime 
we found a VFTH relaxation time above the Kauzmann transition and we derived the 
Adam-Gibbs relation between the relaxation time and the configurational entropy, 
which can be explicitly computed for our model. It is also possible to study the 
dynamics of the system quenched to a temperature below the Kauzmann temperature. 
For such an analysis we refer to . 

We proposed an out of equilibrium thermodynamic formulation and we tested on 
our exactly solvable model whether or not such a generalized approach holds, having 
one extra variable, namely the effective temperature, for the description of the non- 
equilibrium thermodynamics. By effective temperature we mean a thermodynamic 
quantity that would be the temperature of a system at equilibrium visiting with 
the same frequency the same states that the real - out of equilibrium - system at 
temperature T visits on a given time-scale during its dynamics. This kind of parameter 
appears in the thermodynamic functions together with the heat-bath temperature and 
the fields coupled to the system's observables and is coupled to the configurational 
entropy. 

Generally speaking, in order to recast the out of equilibrium dynamics into a 
thermodynamic frame, the history of a system that is far from equilibrium can be 
expressed by more than one effective parameter. This happens when more than one 
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long time-scale is involved in the dynamic evolution of a system. In those cases 
to every time-sector there will correspond an effective temperature Moreover, 
in a given time-sector, the number of effective parameters needed to make such a 
translation into a thermodynamic frame can, in principle, be equal to the number of 
relevant observables considered. In our model, however, for certain dynamic regimes 
determined by the temperature and by the VFTH exponent 7, the effective parameters 
pertaining to processes having the same time scale become equal to each other for large 
times. When the distance from equilibrium is much smaller than the distance 

from the constraint, ^2{t), a single effective temperature alone is enough for a complete 
thermodynamic description of the dominant physic phenomena; this happens when 
7 > 1. 

In many models the concept of effective temperature is used in an attempt to give 
a thermodynamic description of the glassy phase. There are several analytical and 
numerical works defining the effective temperature as the fluctuation dissipation ratio 
(among others |3|, |3^, Q ) and also glassy models in which the effective temperature 
can be measured coupling the system to a slow relaxing thermometer Q . Even 
in granular systems such an approach is followed with some success |36[ . However, 
as far as we know, either no comparison is made between alternative (equally well 
based) definitions of effective temperature or, when such a comparison is carried out, 
very often no coincidence is found H). Nevertheless, in the p-spin model 

the fluctuation-dissipation effective temperature coincides with the thermodynamic 
effective temperature; both are equal to T /x where x is the break point of the Parisi 
function |§ !§. 

This coincidence of effective temperatures for the present model (in those 
parameter regions where it does take place) appears to be connected to a slow enough 
relaxation dynamics. In order to understand what could happen in other models a 
very general analysis should be done to identify what are the minimal requirements 
to produce a thermalization of different degrees of freedom within a single time-scale 
(for a contribution in this direction see, for instance, pO|). 
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